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ABSTRACT: We present the results of a systematic, first-principles study of the spectrum 
and decay constants of mesons for different numbers of color charges N, via lattice com- 
putations. We restrict our attention to states in the non-zero isospin sector, evaluating 
the masses associated with the ground-state and first excitation in the pseudoscalar, vec- 
tor, scalar, and axial vector channels. Our results are based on a new set of simulations 
of four dimensional SU(iV) Yang-Mills theories with the number of colors ranging from 
N = 2 to N = 17; the spectra and the decay constants are computed in the quenched 
approximation (which becomes exact in the 't Hooft limit) using Wilson fermions. After 
discussing the extrapolations to the chiral and large-A limits, we present a comparison of 
our results to some of the numerical computations and analytical predictions available in 
the literature — including, in particular, those from holographic computations. 
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1 Introduction 

The Standard Model of elementary particles describes the strong interaction by Quantum 
Chromodynamics (QCD): a non-Abelian gauge theory characterized by local SU(3) color 
invariance, and with flavors of fermionic matter fields (quarks) in the fundamental rep- 
resentation of the gauge group. Due to asymptotic freedom, the predictions of QCD for 
processes at high energies can reliably be worked out through perturbative calculations. At 
low energies, however, the running coupling becomes large; in this regime many features of 
QCD are determined by non-perturbative phenomena: confinement and chiral symmetry 



breaking. For this reason, the theoretical derivation of low energy QCD properties must 
rely on non-perturbative methods. 

The standard non-perturbative definition of QCD is based on lattice regularization [1], 
which makes the theory mathematically well-defined and amenable to analytical as well 
as to numerical studies. Thanks to theoretical, algorithmic and computer-power progress, 
during the last decade many large-scale dynamical lattice QCD computations have been 
performed at realistic values of the physical parameters, allowing one to numerically obtain 
predictions in energy regimes otherwise inaccessible to a first-principles approach, see, e.g., 
ref. [2] for a recent review. 

A different non-perturbative approach to QCD is based on an expansion in powers of 
1/N, i.e. in powers of the inverse of the number of color charges [3]. When the number of 
colors N is taken to infinity, and the coupling g is sent to zero, keeping the product g 2 N 
fixed, the theory reveals striking mathematical simplifications — see refs. [4, 5] for recent 
reviews. One may study the physical N = 3 case expanding around the 1/N —> limit. 
In particular, in the so-called 't Hooft limit (in which the number of flavors nj is held 
fixed), one finds that the amplitudes for physical processes are determined by a particular 
subset of Feynman diagrams (planar diagrams), the low-energy spectrum consists of stable 
meson and glueball states, and the scattering matrix becomes trivial. Examples of other 
interesting implications of the large- N limit include those discussed, e.g., in refs. [6-10]. 

Another non-perturbative approach to low-energy properties of strongly coupled non- 
Abelian gauge theories is based on the conjectured correspondence between gauge and 
string theories [11-13]. During the last decade, many studies have employed techniques 
based on this correspondence, to build models which reproduce (at least qualitatively or 
semi-quantitatively) the main features of the mesonic spectrum of QCD [14]. Remarkably, 
the large- N limit plays a technically crucial role also in the context of these holographic 
computations: the correspondence relates the strongly coupled regime of a gauge theory 
with an infinite number of colors to the classical gravity limit of a dual string model in 
an anti-de Sitter spacetime, a setup that can be studied with analytical or semi-analytical 
methods. 

In order to understand whether predictions derived from approaches relying on the 
large- N limit can be relevant also for the physical case of QCD with N = 3 colors, it is 
crucial to estimate the quantitative impact of finite- N corrections. For this purpose, recently 
several lattice studies investigated the dependence on N of various quantities of interest — 
including string tensions for sources in different representation, the glueball spectrum and its 
dependence on the 0-term, the topological susceptibility or the finite-temperature equation 
of state — in different SU(iV) Yang-Mills theories [15-34]. These works revealed precocious 
scaling to the large-iV limit: up to modest 0(N~ 2 ) corrections even the SU(3) and SU(2) 
theories appear to be "close" to the large- N limit. Similar results have also been found in 
D = 2 + 1 dimensions [35-43] . 

The purpose of the present paper is to further expand this line of research, by study- 
ing the light mesons in gauge theories based on different SU(iV) gauge groups to explore 
the connections between the different non-perturbative approaches outlined above. More 
precisely, we extract the meson masses from numerical simulations and investigate their 



dependence on the number of colors. We discuss the approach to the large- iV limit and 
compare our results to the theoretical expectations from large- A expansions and holographic 
computations. 

All the results that we present in this paper are obtained from quenched simulations. 
A technical remark about this is in order: in theories including dynamical quarks, virtual 
fermion loops are suppressed by powers of nt/N, so, in general, corrections with respect to 
the large- N limit are expected to be larger than for the Yang-Mills case. However, when A 
is infinite, the theory exactly reduces to a (unitary) quenched model: this justifies studying 
the large- A meson spectrum neglecting dynamical fermions and, in fact, it also provides an 
intuitive explanation of why the quenched approximation performs rather well even for the 
A = 3 case [44]. From the practical point of view, lattice computations of meson masses in 
the large- A limit in the quenched setup offer two major advantages: 

1. neglecting the contribution from virtual fermion loops allows one to bypass the major 
computational overhead associated with the inclusion of the Dirac operator determi- 
nant in the system dynamics; 

2. the convergence to the large- A limit is faster since the leading corrections due to the 
finiteness of the number of colors scale quadratically (rather than linearly) in 1/A, 
enabling reliable extrapolations, even from simulations at rather small values of A. 

In the past, related studies have been reported in refs. [45-49]: one of the goals of 
our present work consists of improving and extending these results, by going to lighter 
quark masses, larger A and larger volumes, and by increasing the statistics, the number of 
interpolators included in the variational basis and of states that we investigate. We also 
aim at clarifying a discrepancy between the results of refs. [45-48], which at large A found 
a value of the vector meson mass close to the one of real- world QCD, and those obtained 
in ref. [49], which, on the contrary, reported a value approximately twice as large. 

In section 2, we define the setup of our lattice computation and, in section 3, we present 
our numerical results. In section 4 we discuss discrepancies with ref. [49]. In section 5 we 
compare our results to analytical predictions and finally, in section 6, we summarize our 
findings and conclude. 

Preliminary results of this work have been presented in ref. [50]. 

2 Setup of the lattice computation 

In this work, we non-perturbatively study theories with SU(A) internal color symmetry 
with A = 2, 3, 4, 5, 6, 7 and 17 color charges, regularized on a finite, isotropic hypercubic 
lattice A of spacing a, in Euclidean spacetime. In the following, we denote the lattice 
hypervolume as L^ x Lt = (Nf x Aj)a 4 . Expectation values of observables are obtained as 
statistical averages over quenched configuration ensembles, using the Wilson discretization 
of the generating functional of the continuum theory detailed below: 
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where dU^(x) denotes the Haar measure for each Uu(x) € SU(iV) link matrix, /j, = 1, . . . , 4 
and S is the "plaquette" discretization of the Yang-Mills action: 
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/3 is related to the bare lattice gauge coupling go via f3 = 2N/g^ and 

£V(x) = U^{x)U u {x + a£)£/J(x + az))C/t(x). (2.3) 

The expectation values of gauge-invariant physical observables O are defined as: 

(0) = ^ [UTl dU ^ 0e ~ s m 

zeAA*=l 

and are numerically estimated via Monte Carlo sampling. The numerical results presented 
in this work are based on sets of configurations generated by code based on the Chroma 
suite [51], using standard local updates [52-57], which we have adapted to work for a generic 
N. In the following, we denote the number of configurations used in our computations (for 
each set of parameters) as n conf . 

Our lattice implementation of quark propagators is based on the Wilson discretization 
of the continuum Dirac operator [1]: 

4 

aM(x,y)[U] = 5 XtV - kJ^ [(1 - i^)U^{x)5 x+a ^ y + (l + 7A«)tfJ(v)$r-aA,i/ » ( 2 - 5 ) 

where the hopping parameter k is related to the bare quark mass m q via: 

(2.6) 

and k c denotes the critical value, corresponding to a massless quark. The additive con- 
stant is given by k^ 1 = 8 + 0(f3~ 1 ) and its non-perturbative determination is discussed in 
subsection 3.1. 

To set the scale, we use the string tension calculations of ref. [58, 59] for N = 2, 3, 4 
and 6, and those of ref. [34] for N = 5 and 7. In each case we choose the coupling 
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such that the (square root of the) string tension in lattice units a\fa ~ 0.2093 is the same 
for each N. 

Using the ad hoc value a = 1 GeV/fm, our lattice spacing corresponds to a ~ 0.093 fm 
or a -1 ~ 2.1 GeV. Strictly speaking, we can only predict ratios of dimensionful quantities. 
In the real world where experiments are performed, rif > 0, N = 3 ^ oo and even the 
string tension is not well defined. This means that any absolute scale setting in physical 
units will be arbitrary and is just meant as a rough guide. 
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Figure 1. A-parameter estimates of eq. (2.8), in units of the square root of the string tension. The 
errors shown are propagated from those of a^/a. 
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4 points 


linear 
quadratic 


208.45 
209.04 


208.16 
208.77 



Table 1. Fit results for f3 at N = 17. 



For the theory with SU(17) gauge group, there are no string tension calculations avail- 
able, so we extracted a /3-value from a fit of the QCD A-parameter in the lattice scheme: 



A « a exp 
with [60, 61] 
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A/^/a was calculated from the data presented in [34, 58, 59] for SU(2 < N < 8) and is 
shown in figure 1. Using the data for N = 6, 7 and 8 and a linear fit in 1/N 2 , we obtained 
(3 = 208.45 for N = 17. Adding further values of N > 4 or using a quadratic fit in 1/A^ 2 
changed this value by less than 0.3% (see table 1). 

Table 2 summarizes the essential technical information of our computations. The N < 7 
results presented in the fits and plots here are obtained from the 24 3 x 48 lattices, corre- 
sponding to L s ~ 2.2 fm. In order to study finite size effects (FSE), we also performed 
additional simulations, both using smaller and (for N = 2 and 3) larger volumes — see the 
discussion in subsection 3.7. In SU(17) we employed a smaller 12 3 x 24 volume. The cor- 
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Table 2. Parameters of the main set of lattice simulations used in this work, j3 denotes the 
gauge action parameter, while k is the hopping parameter appearing in the quark propagator. All 
configurations were separated by 200 combined heatbath and overrelaxation Monte Carlo sweeps 
and found to be effectively statistically independent. For orientation we also include the bare 
't Hooft parameter. 



responding extent L s = 12a ~ 2.512 /y/a is well above the inverse critical temperature [34] 
T c 7 l « 1.681/V^ < L c . 

The K-values were selected to keep one set of six pion masses approximately constant 
across the different SU(iV) theories. To achieve this, we combined the results reported in 
ref. [46] for the groups studied therein with initial estimates for the groups that had not 
been studied before. We vary the "pion" mass down to m^/y/a ~ 0.5 for groups with 
N > 5, and to m^/y/a ~ 0.75 for N < 5. We also simulated a smaller quark mass for 
SU(A r < 5) but found significant numbers of "exceptional configurations" [62] (up to 15 % 
of the total); we leave these data out of this work. For N = 5, at the lowest quark mass, 
only two exceptional configurations were encountered that we removed from the analysis. 



2.1 Smearing and operators used 

The meson spectrum can be extracted from zero-momentum correlators of interpolating 
operators of the form: 

O k (x,t) =tf(x,t)T k if>(x,t), (2.10) 



where different choices of T^ correspond to different physical states. We used degenerate 
quark masses for the quark fields u and d, which allowed us to study the spectrum of the 
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Table 3. List of the studied channels and their bilinear operators used in the correlation functions. 

particles n, p, ao,ai, b\ (see table 3). A generic meson correlator can then be computed as: 
Crr(t) = ]T(O(x,t)O(0)) = -^Tr(rG(x,t)r 75 G(x,t)75) , (2.11) 



where the propagator G(x, t) is obtained by inverting the Dirac operator with the stabilized 
biconjugate gradient (BiCGStab) algorithm with even/odd preconditioning. The trace 
indicates a sum over spin and color. In the rare cases, when the BiCGStab routine failed 
to converge (which only occurred at the lowest quark masses), we reverted to the standard 
Conjugate Gradient (CG) algorithm. 

In our correlators, we used both point- like and extended sources and sinks, employing 
several steps of Wuppertal smearing [63] which iteratively modifies a fermion fields as: 
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where n denotes the number of iterations, uj is the smearing parameter (we used u = 0.25) 
and U' is the gauge field smeared by 10 iterations of the spatial APE smearing routine [64]: 



U-(x) = Proj su(7V) a Ui(x) + ^ Uj(x)Ui(x + aj)U^(x + 
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with smearing parameter a = 2.5. 

2.2 Variational method 

We extracted the ground state and the first excited level using the variational analysis 
discussed in ref. [65-67]. For each channel, we computed the cross-correlation matrix 
Cij(t) = {Oi(t)Oj(0)), where % and j correspond to the number of iterations (0, 20, 80 
or 180 steps) of Wuppertal smearing, for the sources and the sinks. 
Then we solved the generalized eigenvalue problem: 



C(t)v a = A a (t)C(t )v c 



(2.14) 



and extracted the mass m performing hyperbolic-cosine fits of the largest and second largest 
eigenvalues, A and A 1 , in [t min ,tmax] ranges (t max < aN t /2): 



\(t) = A (e~ mt + e -m{aN t -t)\ _ 



(2.15) 



All statistical uncertainties were estimated using a jackknife procedure. With four different 
operators Oi in many cases we were able to extract the first three states, however we regard 
the second excited states as unreliable at the present statistics. Using a subset of three 
operators out of the four mentioned above leads to compatible mass values within errors. 

Varying to in the range [0, 2a] gives compatible results, so we used to = a. To select 
the best fit ranges, for each particle we first studied the effective mass defined as: 

'X(t + a) + \{t-a 
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2A(t) 

We determined £ m j n as the Euclidean time separation at which m e s reaches a plateau (within 
statistical uncertainties), so that the contribution from higher states is negligible, while t max 
is the value where m e ff becomes too noisy for stable fits. The signals become more precise 
at larger N and lower K-values. Typically we fit A(i) in the range [5a, aiVj/2] for the ground 
states and in the range [5a, 10a] for the first excited states, adjusting those ranges (by one 
or two lattice spacings) on a case-by-case basis. Fitting to eq. (2.15) with this procedure 
leads to reduced x 2 - varues which are well below one. 

The cost of inverting the propagator increases at lower quark masses, with the signal 
becoming noisier at the same time. For this reason we focused only on the ground states 
for the lowest two quark masses of each SU(iV) group and instead of using the variational 
method we computed the two point functions using only 80 steps of smearing for the 
sources /sinks. We then applied the same analysis for \(t) directly to the correlator C(t). 

3 Numerical results 

3.1 The PCAC mass and the critical hopping parameter k c 

Our non-pert urbative determination of k c is based on the partially conserved axial current 
(PCAC) relation, 

X>|V^(F,t)|7r> = 2m PCAC (t) X>|P(r,i)|7r), (3.1) 

r r 

where A fl (r, t) = u(f, t)j^ 5 d(r, t), P{r, t) = u(f, t)jnd(f, t) and m PCAC = lim^oo m PCAC (t). 
On the lattice we compute m PCAC as: 

„™ f+\ — ^7075,75 (^ + a ) ~ ^7075 ,75 \t ~ a ) /o o\ 

where the pion sources are smeared, and fit this to a constant. This quantity, which is not 
affected by chiral logarithms, is related to the vector quark mass am q defined in eq. (2.6) 
by renormalization constants: 

WPCAC = ry ry "V (3.3) 

Taking into account the leading lattice corrections {l + bx^m) to renormalization constants 
Zxi we can fit our lattice results to the expression: 

Z P . 1/1 1 \ , A , 

a TOpcac = ry v (1 + Aa m PCAC + ...)- , (3.4) 






from which we extract k c , Zp/{ZaZs) and A (see table 8 of appendix B ). The fits are 
plotted in figure 2. 

The (unrenormalized) PCAC mass can be determined very precisely (see tables 10-16), 
so we will expand every meson mass as a function of this variable. 

We expect the parameters A, Zp/(ZaZs) and k c to have 0(1/N 2 ) corrections, hence 
we fit them to 

"1 + ^2, (3-5) 

as shown in figure 2. With this analysis, we obtain good fits for A and Zp/(ZaZs), with 
values of the reduced \ 2 close to 1, while for k c , although qualitatively the behavior looks 
very promising, we get a \ 2 P er d.o.f. of 300, indicating that uncertainties in our data are 
underestimated. In fact k c can be considered a function of (3 only, which was chosen to 
match the string tension among the different groups. This process introduces a systematic 
error which propagates to k c . A qualitative way to estimate this propagation is to consider 
finite differences between the values of k c ,/3 and (i\fo from refs. [34, 45, 58, 59] and to 

compute 

An c A/3 



5Kc = "AflATr^' 5 ^- ( 3 - 6 ) 
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The r.h.s factors of the equation above are listed in table 9, together with 5k c for the groups 
available in refs. [45, 46]. These systematic uncertainties due to the matching of j3 are ten 
to twenty times larger than the statistical ones and, since they are approximately constant 
across the N- values, we used the same errors for the remaining groups. Taking this into 
account, the reduced x 2 -value of the k c fit becomes 1.6. 
The 1/N 2 fit results are: 

Z P , x 0.699(45) . , 

- 0.8291(20) -^-L, (3.7) 



Z A Z S v ' N 2 

A = 0.390(13) + ^M, (3.8) 

, w , 0.028242(68) (394) . , 

k c = 0.1598555(33)(447) -^-^ -, (3.9) 

where in the K c -case the second error is the systematic one, due to the slight mismatch 
in the string tension, detailed above. We find the ratio Zp/(ZaZs) to vary between 0.68 
(N = 2) and 0.83 (N = 17), with the SU(3)-value 0.75, which is consistent with the 
non-perturbative result 0.82(11) [68] obtained at /3 = 6.0, close to our value /3 = 6.0175. 

3.2 The pion mass 

The pion masses are shown in figure 3 and presented in tables 10-16 of appendix B as 



dimensionless ratios, dividing them by the square root of the string tension y/a. Quenched 
chiral perturbation theory [69] predicts 

{am n ) 2 = A(am q )^+i + ..., (3.10) 

where the exponent 5 is due to the presence of chiral logarithms. The theory predicts 5 to 
be positive, O(10 _1 ) for SU(3), and suppressed as 1/N at large N. However, it is known 
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Figure 2. Fit of the PCAC mass to eq. (3.4) (left), together with the 1/N 2 fit of the parameters 
A, Zp (ZaZs)~ and k c (right). We plot systematic errors for k c . 



that data in the region where m^/y/a > 1 are not very sensitive to chiral logarithms [70], 
and fitting larger pion masses according to eq. (3.10) would lead to values of 8 with even the 
wrong sign. For this reason, we included a subleading term of the quark mass expansion, 
performing fits according to: 
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(3.11) 



In these fits, the 5 exponent is, essentially, determined by the lowest pion masses — which, 
unfortunately, are the points with the largest uncertainties. This leads to rather large 
relative errors for 5. Nevertheless, we found clear evidence that S gets smaller when N is 
increasing. Within our precision limits, 5 is found to be consistent with zero for all N > 6. 
In fact, for larger N one can omit 5 completely from the formula and still obtain a good 
fit. Conversely, for N < 3 we find ©(lO^ 1 ) values of 5, where for SU(3) we get results 
consistent with those reported in ref. [71], in which Wilson and clover actions were used, 
and also consistent with table 3 of ref. [70], where different values of 5 were calculated, 
using different actions. This suggests a small 1/N coefficient and thus we include a higher 
order term into the 1/N fit. 

The N dependence of 5 and the expansion of A and B in powers of 1/N 2 (see figure 3) 
give a x 2 P er degree of freedom close to 2 and read: 
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In order to assess the systematics on the exponent S, we performed a combined fit 
(m PCAC ,iV) of our data to eq. (3.11), using all the iV-values at once and excluding the 
two highest masses for each group. Since the Wilson action explicitly breaks the chiral 
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Figure 3. Fit of the squared pion mass, in units of the string tension, to eq. (3.11) (left panel), N 
dependence of the fit parameters eq. (3.12)-(3.14) (right). 



symmetry, we include in the fit also a constant term which however is found to be consistent 
with zero. The resulting curve has x 2 /d.o.i. = 1.6 and reads, 

,2 



a 



0.0015(36)+ 11.67(15) 



8.1(5.4) 

N 2 



m F 



a 



i 

1+6 



+ 



+ 2.95(42) 



1(15) \ m 



""PCAC 



a 



0.093(27) 1.00(52) 



N 



AT3 



(3.15) 
(3.16) 



Note that the m 2 CAC term is now less well determined, due to the exclusion of high 
mass points. However, the 5-parametrisation is consistent with the one of eq. (3.14). 

In figure 4 we plot for each group the pion mass according to eq. (3.15) divided by 
the PCAC mass, in order to emphasise the deviations from a linear behaviour, due to the 
exponent 5. 

Below, we expand all the remaining meson masses as functions of m PC Ac- These can 
easily be translated into dependencies on m 2 through eq. (3.11) above. 

3.3 The p mass 

Quenched chiral perturbation theory predicts a dependence of m p on the square root of the 
quark mass m q [72], in contrast to the unquenched theory, where the leading behavior is 
linear in m q . Thus the expansion of m p takes the form: 

rn p = m pfi + C 1 / 2 m 1 q /2 + C\m q + C 3 / 2 m 3 q /2 + . . . , (3.17) 

where the C]/2 coefficient is expected to be negative and to vanish as 1/iV in the large- N 
limit [72], restoring the linear behavior. 

The masses of the p states are listed in tables 10 to 16 and plotted in figure 5 against 
the PCAC mass. For each group we fit: 



m 



A + B 



m 



PCAC 



a 



1/2 



+ c ,mpcAc +jD 



a 



m 



PCAC 



a 



3/2 



(3.18) 
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Figure 4. Combined fit of rri%./ {m VCAC *Ja) . The solid curves are calculated using eq. (3.15), for 
each N. 

and then expand A, C and D in powers of 1/N 2 , while fitting B linearly in 1/N . Note that 
D should vanish like 1/N 2 in the large- N limit. The parameter B is found to be compatible 
with zero for N > 5 while the parameter D has unexpectedly a finite value: this effect is 
due to contamination from higher orders in the quark mass, as we show below. 
The large- A r expansion of the parameters reads: 

1/2 



TO, 



a 



1.504(5!) + H3P2) _ ™M 



+ 3.08(53) + 



TV 2 
16.8( 



TOr 



a 



8(8.2) \ / topcacA 

iV 2 U? 



-0.84(31) 



9.4(4. 



3/2 



8)\ f m PCAC \ ' 



19) 



A^ 2 J V V^ / 
To address the question of the non- vanishing large- N value of D, we interpolated the 



data with an alternative fit (figure 6) of the form: 

0.51(11) 



to,, 



a 



1.5382(65) + 



vY 2 



+ 2.970(34) 



3.39(55) \ 
iV 2 J 



m 



PCAC 



-0.706(43) + 



3.00(68) \ / topcac V 



V v " / 



(3.20) 



Although the 1/N counting of eq. (3.20) is not consistent, because of the missing square root 
term, one should notice that eq. (3.19) and eq. (3.20) share the same N — > oo behaviour. In 
fact the N = oo coefficients of eq. (3.19) agree within the errors with those of eq. (3.20). In 
particular, we stress that the non-vanishing term 0.84(31) (m PC Ac/V^) 1S consistent with 
0.706(43) {m-p CAC /^/a) , for the range of m PCA c studied. In principle one should introduce 
a quadratic term into eq. (3.19), in practice however, such a fit to six points with five free 
parameters becomes unstable and does not allow us to study the 1/N behaviour of the 
coefficients. 
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Figure 6. Quadratic fit of the p mass in units of the square root of the string tension, according 
toeq. (3.20). 

A possible solution is to fit the data for all groups at once using a combined fit, where 
we fix the (m PCAC , N) functional form up to the second order term in the quark mass 
(figure 7). With this approach we obtain: 



rrir 



a 



+ 



1.5395(83) + 
27(11 



0.92(21) \ 0.06(14) 



N 2 



N 



w T 



a 



1/2 



+ 2.994(44) 



13.9(48) 

N 2 



m T 



m 



PCAC 



a 



3/2 



+ -0.739(50) 



15.1(73) 

N 2 



m 



PCAC 



a 



(3.21) 



with a x /d.o.f. = 2. The drawback of this approach is that data at small (m q ,N) have 

1/2 

less weight in the combined fit, leading to a smaller coefficient for the m PCAC term. 

In order to compare our results with the holographic prediction (see section 5.1) it 
is useful to study m p as a function of the pion mass, paying particular attention to the 
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Figure 7. Combined fit of the p mass in units of the square root of the string tension. The solid 
curves are calculated from eq. (3.21) for the different N. 

linear term of m p {m^). Since the results of eqs. (3.19)— (3.21) tend to agree in the large- N 
limit and considering that our extrapolation relies on small N, which might be affected by 
quenched deviations, we quote eq. (3.19) as our best phenomenological parametrization of 
the data. To avoid further propagation of errors, we fit for each JV the p data directly to 



the 7T masses using the fit form of eq. (3.19), with the substitution m PCAC 
extrapolate the fitted slope to N — > oo, obtaining: 



m 



t. Then we 



mptrriir) 



1 + 0.360(64) 



m. 



+ .... 



(3.22) 



m p (0) ' \m p (0) 

As a consistency check, we can omit the data corresponding to the largest quark masses 
and repeat the analysis using only the smallest three masses and a simple linear fit. This 
leads to: 



771^(771*.) 



1 + 0.317(2) 



771-7. 



(3.23) 



m p (0) ' \m p (0) 

where the smaller error may be unreliable because the fit has two free parameters for only 
three data points. We quote eq. (3.22) as our final result. As already pointed out in ref. [46], 
the results for the slope are close to the holographic prediction obtained in the context of 
the model presented in ref. [73] (see ref. [14] for a detailed discussion). 

3.4 The scalar particle 

The analysis of the scalar mesons ao requires special attention. In the quenched theory, in 
which 7/ also becomes a Goldstone boson violating unitarity, the scalar correlator shows a 
long-range negative contribution, in addition to the standard short-range exponential decay. 
In ref. [74] it was shown that this effect, which is dominant and clearly visible only at the 
lowest quark masses, is caused by loop diagrams corresponding to an intermediate rf — n 
state, which is light and has negative norm in the quenched approximation. Our approach 
consists in fitting the ao two-point functions as: 

C(t) = C e- ma o* _ Cie-i", (3.24) 
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where the "unphysical" quantities C\ and /x are fixed using the values at large t only. This 
approach works very well for N > 5, where the noise in the central region of the correlator is 
smaller. For similar smearing/normalization we expect the amplitudes C, to be proportional 
to N, with 1/N 2 corrections, 

Ci^NLi + ^ + .-.Y (3.25) 

In particular, we expect the ratio C\/Cq to be zero at N = oo, i.e. a\ = and C\ oc iV -1 . 
Indeed we find no evidence of negative contributions at N = 17 for the values of the quark 
masses studied. Moreover, the negative contributions get smaller at higher masses, meaning 
that C\ is at least suppressed like mp^, AC (and vanishes for m ao > 2m n ). (Un)fortunately 
this behaviour is visible only at the very lowest quark mass for each of the groups studied, 
shown in figure 8 (left), so we cannot analyse in detail the amplitude C\ as a function of 
m PCAc- The best strategy to give a qualitative estimate of C\ is to employ a combined 
(m P cAC) N) fit, where we use data points from all group (figure 8 - right). This leads to the 
preliminary estimate 



0.0269(53)—^ , (3.26) 

with a reduced x 2 of 1.4. Note that for SU(2) (not shown in the figure) we were unable to 
obtain meaningful results. 

The ao masses, calculated according to eq. (3.24), are listed in tables 10-16 and their 
1/N 2 expansion (plotted in figure 9) reads, 

""<» - ^2.402(34) 4- i^>) + (2.721(53) - «^fi) I2S2. (3^27) 



a 



3.5 The remaining mesons 

The ground state energies for the remaining states are listed in tables 10-16. The corre- 
sponding fits are shown in figures 15-16 of appendix B; in these cases, our fits do not include 
quadratic terms, because of the larger uncertainties of the data. We fit the remaining states 
to the form, 

mx _ ( A , A xA,f„ , B x , 2 \ Wpcac , , 

7^ - \ Ax ' 1 + ~w j + 1^- 1 + -n 2 - -^rT (3 ' 28) 



and list the results for Ax,i,Ax,2)Bx,i and Bx,2 in table 24. 

As discussed above, the variational method allows us to extract the first excited states in 
each channel, although the results for these masses are much noisier than the corresponding 
ground states. We list the masses in tables 17-23 and plot in figures 17 to 21 their quark 
mass dependences and the 1/N 2 dependences of the respective fit parameters. The excited 
states are also fitted to eq. (3.28) and the results are listed in table 24. 
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A and B (right). 



a = A+B-mpcAc/y/v (left) and 1/N 2 fit of the parameters 



3.6 Decay constants 

We define the lattice meson decay constants 1 F^ and /' at as 

(0|A 4 |tt) = V2m n Fl a \ 



\V k \p\) = m p ff % e k {p,X), 



(3.29) 
(3.30) 



where A^ = 147475^ and Vk = ujkd are the non-singlet axial and vector currents, respec- 
tively, while ek(p, A) denotes a polarization vector. 



1 Throughout the paper we use the notation F n — / 7r /v / 2 ~ 92 MeV to indicate the pion decay constant 
and Fl at to address the same quantity calculated on the lattice up to the renormalization factor Za, i.e. 
F, = ZaF 1 ^. 
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We fit the pion and the axial correlators for large t as: 

(A 4 (t)n a (0)) ~ ^L( |A4|7r)(7r|7rt |0)e- tm - = C Ai e~ tm " (3.31) 

(7Ta(t)7r a (0)) x ^-(0|7r Q |7r)(7r|7rt |0)e- tm ^ = C^~ tm - (3.32) 

(where 7r Q is one of the four differently smeared pion interpolators) and then compute F^ 
as: 

In order to better compare large- N to N = 3 results, we choose to normalize the decay 
constants as: 

For the SU(3) theory, we find results which are consistent with ref. [75], while for larger 
gauge groups our measurements show that the expected \/N scaling behavior is well satis- 
fied, see figure 10 and tables 25-31. 

The 1/N 2 expansion of the (rescaled) pion decay constant reads: 

Fl at ( , x 0.121(56) \ / , , 0.29(30) \ m PCAC 

^0.2619(37) - -^ r A ) + (0.506(24) - -^-i J -^f 



f « = f *\Nf> U = fp\hf- ( 3 - 34 ) 



(-0,20(3!) + «)^ 



+ (-0.320(31)+ ;;; J ) -^. (3.35) 

In the case of /' at , we use a similar approach — the main difference being due to presence of 
the polarization vector e^(p, A). This satisfies the relation 

^2e^(p,X)e u (p,X) =gp,v- -A^-, (3.36) 

A P 

so that for zero momentum and for a fixed direction \i = v the above expression becomes 
one, and the computation is identical to the previous one (up to a \/2 factor). To improve 
the statistical precision, we averaged the results over the three spatial directions. We plot 
the p decay constants for different JV in figure 11; the 1/N 2 fits can be summarized as: 

% = (0.8173(70) + ™) + (0.582(44) - ™) "^ (™V 

2 



+ (-0.407(57) + i^2)=^c. 



The decay constants computed on the lattice are related to the ones in the continuum by 
the renormalization constants Za and Zy, 

F w = Z A F™, f p = Z v f l p at . (3.38) 

Za and Zy have been determined non perturbatively only for N = 3 [68, 76] while two-loop 
perturbative results are known [77, 78] to converge slowly. In appendix A we detail our 
strategy to obtain estimates of Za and Zy for N ^ 3 and of their systematics. The results 
for F n and f p are listed in tables 25-31. 
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Figure 10. Fit of the rescaled pion decay constant to: F^/y/a = A+5-to PC ac/v / ^ : +C' to pcac/ iT 
and 1/N 2 fit of the resulting fit parameters A, B and C ' . 





0.5 




1 1.5 1.75 2 


2.25 


2.5 


2.75 




- 1 ' 


1 ' 1 ' 


1 


1 ' 1- 

T 




— nx 


' 








^^ggsS^--^^ 


ZZ-^£ 




0.95 


— 


T 




^ — -~ 




- 








'f\^^ """ 






j - 




i 


--' 






_ 


0.85 




SU(2) 
SU(3) 




■?■ 








SU(4) 




0.8 


z 








SU(5) 
SU(6) 
SU(7) 
SU(17) 


- 


0.73 


' 




1 




-- SUM 


I 






, 










1 ' 1 


- 


-0.2 








- 


-0.4 








- - 


-0.fi 


i 




1 


- 


0.6 


T T 




1 


- 


0.4 


1 






■ - 


0.2 






1,1- 


—^ 


0.88 


T 




1 


- 


0.84 


1 


T 


T 


- - 


o.x 


: h± 




i , i ■ 


- 



Figure 11. Fit of the rescaled p decay constant to: fp/yfcr = A + B ■ m PC Ac/ \l& + C ■ rrip CAC /a 
(left), together with the 1/N 2 fit of the resulting parameters A, B and C (right). 

3.7 Finite volume effects 

Due to the Eguchi-Kawai volume reduction [79], finite size effects (FSE) are expected to be 
zero at infinite N, as long as all lattice dimensions (in physical units) are kept larger than 
a critical length L c [80], so that center symmetry is not spontaneously broken. 

At finite N, FSE become larger for smaller quark masses and for smaller N [81, 82]: 



m^L) = m^oo) [1 + B exp(—m n (oo)L))], 



(3.39) 



where the parameter B vanishes in the large-i¥ limit. For N = 2 and 3, we carried out 
simulations at three volumes and fitted the pion masses to eq. (3.39), obtaining the results 
displayed in figure 12. As one can see from these plots, FSE are drastically reduced going 
from 2 to 3 colors, where the data can already be fitted to a constant. For larger A^-values 
we carried out simulations at two volumes only, unsurprisingly, without any evidence of 
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Figure 12. Fit of m„(L) /m w (oo) according to eq. (3.39) for the SU(2) (left) and SU(3) theory 
(right). 

FSE. While one may also carry out a similar analysis for the other particles, their finite 
size corrections are expected to be smaller than for the pion, and thus negligible within our 
statistical uncertainties. 

3.8 Finite cut-off effects 

Before comparing our results to theoretical continuum predictions, it is important to discuss 
the potential impact of discretization effects on our lattice results. Strictly speaking, since 
our simulations were performed only at one value of the lattice spacing, it is not possible 
to perform a continuum extrapolation. However, following the discussion in ref. [46] (a 
study carried out at values of a very close to ours), one can nevertheless get an estimate 
of the systematic uncertainty induced by the finite lattice cut-off, by comparing the results 
with those obtained at a different, coarser, lattice spacing in ref. [45]. In particular, in 
ref. [46] it was shown that, in spite of the 60% difference between the two lattice spacings 
used in the two studies, the p meson masses obtained are very close to each other (up to 
differences of the order of 5%). Since our simulations (like those reported in analogous 
previous studies [45-48]) use the standard, unimproved, Wilson discretization for the Dirac 
operator, the leading lattice artifacts are expected to scale like a (those in the gluonic 
sector are proportional to a 2 , and, hence, subleading). For the p meson mass, the analysis 
carried out in ref. [46] showed that, by extrapolating their results and those from ref. [45] 
either linearly or quadratically in a, one obtains estimates of the continuum-limit value 
which are about 2% to 5% off from the result obtained at a ~ 0.09 fm. In the absence 
of analogous data for the other states investigated in the present study, it is reasonable to 
take these numbers as a rough, order-of-magnitude estimate of the systematic finite cut-off 
errors affecting our results. 

As a final remark, we observe that while taking the continuum limit can lead to a slight 
quantitative change in our results for physical quantities, it should not dramatically affect 
their (mild) dependence on N. In principle, interchanging the order of the continuum and 
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large- N limits of a lattice model may involve some subtleties (see ref . [83] for a discussion) . 
However, the results of previous studies indicate that, for lattice gauge theories and for the 
observables discussed here, the two limits commute [4, 5]. 

3.9 Large-iV spectrum 

To conclude this section, we display in table 4 the results for the meson spectrum extrapo- 
lated to infinite N at different quark masses. The results are listed in units of yfa and in 
units of the (normalized) pion decay constant in the chiral limit, 

Foo = \J^F«{rn q = 0) , (3.40) 

which should be particularly useful for chiral perturbation theory (xPT) applications [84- 
86] (see section 5.2). 

Of phenomenological interest is not only the spectrum at m q = (figure 13) but are 
also the spectra at m q = m u d and at m q = m s (figure 14) where m U( i and m s denote the 
physical (isospin- aver aged) light quark and strange quark masses, respectively. We fix the 
former imposing at N = oo the values [87]: 

Foo = 85.9 MeV, (3.41) 

m^rriud) = 138 MeV. (3.42) 

Foo fixes the lattice spacing to a = 0.10 fm and the string tension to a = (395 MeV) 2 . This 
value leads to a better agreement between the large- N meson spectrum and the experimental 
masses than the previously used estimate a = (444 MeV) 2 (see figure 14-left). We find the 
ratio 

F 4 m ^ = 1.020(20) , (3.43) 

F w (0) 

where the renormalization constants cancel, to be compared to the value F n -(m u d)/Fr(0) = 

1.073(15) for N = 3 QCD with sea quarks [87]. In the large- N limit the mass-dependence 

of F n appears to be reduced. 

The strange quark mass is obtained by fixing at N = oo the mass of a (fictitious) 

strange-antistrange pion to the value 

m n (m s ) = (m 2 K± + m 2 R0 -m 2 ± ) 1/2 « 686.9 MeV. (3.44) 

The large- iV meson spectrum at m q = m s is shown in figure 14 (right). 

Note that our way of fixing m uc i and m s is arbitrary and different choices of input 
observables of real N = 3 QCD may give values that differ by 0{\/N) corrections. In fact 
the procedure of adjusting the quark masses described here differs from that of ref. [50]. 

Setting the scale from Foo = 85.9 MeV and fixing the quark masses as detailed above, 
we find the values m p = 619(9) MeV at m q = m u d and 876(13) MeV at m q = m s , to 
be compared to m p = 775 MeV and m^ = 1019 MeV in QCD. For the scalar mesons 
we find m ao = 959(18) MeV and 1208(22) MeV for light and strange quarks, respectively. 
Note that the errors stated do not include an overall scale uncertainty of 8 %, due to the 
renormalization of the pion decay constant. 
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Figure 13. The meson spectrum for different N in the chiral limit. The masses and the decay 
constants are given in units of the square root of the string tension for each group SU(iV) and the 
extrapolated N — >■ oo values are shown as horizontal bands. 




4.6 1 







Ground stt 

Lac ill.' J ski 



Figure 14. Left: the large- N meson spectrum at the m q — m u ^ scale compared to the experimental 
2V = 3 masses for the two values of the string tension: a — (395 MeV) 2 , obtained from imposing 
Foo — 85.9 MeV, and the ad hoc value a — (444 MeV) 2 . The pion mass is taken as input to set 
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m u d- Right: the large-V meson spectrum at the m„ = m s scale. 



4 Comparison with other lattice studies 

The mass of the p meson has generated some controversies in the literature: refs. [45-48] find 
a value that is close to the one measured in SU(3), while ref. [49] finds a mass approximately 
twice as large. In this work, we confirm the results of refs. [45-48]. Nevertheless, the 
discrepancy with ref. [49] needs to be taken very seriously and properly addressed, since 
the results of ref. [49] question not only the value of the p mass, but (indirectly) the whole 
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2.428(31) 


3.059(33) 


11.04(21) 


11.17(21) 


14.07(25) 


a\ 


1++ 


2.860(21) 


2.883(21) 


3.414(23) 


13.16(21) 


13.26(21) 


15.71(24) 


h 


1+- 


2.901(23) 


2.924(23) 


3.452(24) 


13.35(21) 


13.45(22) 


15.88(25) 


vr* 


o-+ 


3.392(57) 


3.413(57) 


3.887(60) 


15.61(34) 


15.70(34) 


17.88(37) 


P* 


1— 


3.696(54) 


3.714(54) 


4.127(56) 


17.00(34) 


17.08(34) 


18.99(37) 


ao 


0++ 


4.356(65) 


4.375(65) 


4.816(69) 


20.04(41) 


20.13(41) 


22.16(44) 


at 


1++ 


4.587(75) 


4.605(75) 


5.012(81) 


21.10(46) 


21.18(46) 


23.06(49) 


b\ 


1+- 


4.609(99) 


4.628(99) 


5.06(11) 


21.20(54) 


21.29(55) 


23.29(58) 


F„ 


- 


0.2174(30) 


0.2216(30) 


0.3035(58) 


1 


1.020(20) 


1.396(33) 


u 


- 


0.5721(49) 


0.5762(49) 


0.6516(93) 


2.632(43) 


2.651(44) 


2.998(60) 



Table 4. The N = oo meson spectrum and decay constants in units of the square root of the 
string tension y/a and in units of the (normalized) chiral pion decay constant F^ = F w (0) \JS/N 
in the chiral limit, for three different values of the quark mass which are set using the pion masses 
as inputs. We also display i*V = F^ y/3/N and f p — f p y/3/N in the last rows of the table. The 
results have not been extrapolated to the continuum limit. To account for this, based on [46], a 
systematic error of 5% should be associated to all values. Regarding the pion decay constant, due 
to the lack of non-perturbative renormalization at N = oo, an error of 8% should be associated to 
the values in units of i 7 ^ and to the decay constants in units of the string tension. 

meson spectrum obtained in this work. The issue is not straightforward to settle, since the 
method used to address the problem in ref. [49] is different from the one in refs. [45-48]. 
Another related issue concerns F n (N = oo) [88-91] that was used in [49] to set the physical 
scale. Our result, indicated in table 4, is 



F^/xfo = 0.2174(30) 



(4.1) 



while from [88] one infers 



^00/^ = 0.277(12) 



(4.2) 



where the known value [34] T c /y/a = 0.5949(17) has been used to convert the original 
results (provided in units of the deconfinement critical temperature T c ) into units of \fa. 
Note that in both estimates we have considered only the statistical error. When also the 
systematic errors related to the finite lattice spacing and (for the calculation described 
in this work) to the renormalization constant are taken into account, the results agree 
within two standard deviations, showing an acceptable level of consistency. Instead, the 
discrepancy of the results for the mass of the p meson demands a thorough inspection of the 
numerical techniques and of the calculations that lead to the conflicting numerical values. 
Refs. [45-48] use the conventional techniques we described in section 2 to extract the 
mass of the p from the long-distance behavior of correlators of fermion bilinears carrying 
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the quantum numbers of the state of interest, at relatively low N (up to 7); this method 
has been well tested in standard QCD calculations. On the other hand, the authors of [49] 
extract the mass at larger N (N = 17 and N = 19), but using smaller lattice sizes, from 
momentum space correlators which are obtained injecting allowed discrete momenta in the 
lattice configurations. Although the quenched momentum technique used in [49] is less well 
tested, one cannot exclude the possibility that it gives better results at large N. In addition, 
the fermion discretization is different: Wilson fermions in [45-48] and in the present work, 
and overlap fermions in [49]. 

Taking the technical differences into account, in principle the discrepancy may have 
several interpretations: 

• it may be due to large artifacts in either of the lattice computations (which are 
performed at finite quark masses and lattice spacings, and use different discretizations 
for the Dirac operator — characterized, in particular, by different chiral properties); 

• the quenched momentum computation in refs. [49] is based on an expected cancellation 
of finite-volume effects in the large- N limit. If such cancellation is not complete, the 
numerical results could have been contaminated by the associated systematic effect. 2 
A (possibly related) issue could be the mixing with some excited state(s), given that 
ref. [49] extracts the meson spectrum by computing quark propagators at a small 
number of low momenta; 

• large- N gauge theories have a complicated phase structure [80], which might create 
metastabilities. It is hence possible that the two calculations are in different phases, 
with (possibly) only one being in the phase connected with the infinite volume limit; 

• the discrepancy may actually indicate a "physical" effect: the results reported in 
ref. [49] are obtained at much larger values of N than those investigated in refs. [45- 
48], and it may happen that, for a moderate number of colors, the finite- N correction 
coefficients conspire to hide the deviations from the large- N limit, which might become 
visible only for much larger values of N. This interpretation, however, is at odds with 
the very smooth ./V-dependence of all other observables that have been investigated 
so far. 

Since it is of paramount importance to resolve the root causes of the discrepancy, we 
performed a dedicated SU(17) calculation at values of the parameters that are very similar to 
those of ref. [49] (/? = 208.45 vs. /3 = 208.08), to enable a direct comparison. Interestingly, 
setting the scale by the string tension, it turns out that refs. [46, 48] and the present study 
use a similar lattice spacing as [49]. From this data set, we obtained values of the p mass 
that are compatible with our results at lower N. Finally, we have monitored the average 
local value of the Polyakov loop and of the plaquette in the four directions, to check that 
the whole lattice is in the phase relevant for the continuum limit of the large- N theory. 
One might object that since fluctuations are reduced at large N, extracting masses with 

2 Although, as we have already discussed, we do find some level of consistency for Foo, it is possible that 
the technique has different systematics depending on the quantum numbers of the meson. 
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connected correlators might not be the correct strategy. This issue has been addressed 
in [46] , where it has been shown that for mesons the relative precision of the mass in fact 
improves as N increases. 

As a consequence of this analysis, where we find N = 17 results to be in good agreement 
with N < 7, the possible sources of the discrepancy can be narrowed down to the different 
fermion discretization, the different method for extracting the mass and the influence of the 
bulk phase on the numerical results. Concerning this last point, we remark that while in our 
calculation we have kept fixed a physical scale (provided in our case by the string tension), 
in [49, 88] /3 is varied across the gauge groups using the large-iV perturbative /3-function. 
This might underestimate the value of /3 at which the bulk is present at large- iV (one can 
compare the f3 used in [88] with the investigation of the bulk phase reported in [58] to 
realise that there could be an influence of this unwanted phase on some of the results of 
[88]). This could be the source of the small difference in F^. In fact, it was noticed already 
in [49] that b = (3/(2N 2 ) = 0.35, which was crucial for obtaining the value of F^ in [88], 
gives a variation of the mass of the p of around 30% with respect to the value obtained at 
b = 0.36 (/3 = 208.08), and this variation was ascribed to bulk effects. It is then reasonable 
to assume that the same effects might have influenced the extraction of F^ in [88] . 

To sum up, while the techniques described in this work are well established, with 
known systematic effects (which we have discussed in detail), those of [49, 88] are new 
and further work is needed before the reliability of the corresponding results can be fully 
assessed. Should future investigations confirm them, clearing the doubts raised in this 
paper, the implications will be instructive and a new intriguing picture of the relationship 
between QCD and the large-iV limit might emerge. In the absence of such an analysis 
and supported by the numerical results for N = 11 obtained in this work, the numerical 
evidence that we have presented shows a large- N meson spectrum and decay constants that 
are close to the SU(3) ones, with the N = 3 values described by modest 1/N 2 corrections 
to the observables at N = oo. 

5 Comparison with analytical predictions 

In this section, we relate our numerical results to analytical predictions. We mainly focus 
on calculations derived in the context of gauge/gravity models (subsection 5.1), and briefly 
mention some of those in chiral perturbation theory (subsection 5.2). 

5.1 Holographic models 

Part of the motivation for studying the meson spectrum in the large- iV limit comes from 
the Maldacena conjecture [11-13], namely from the expectation that gauge theories admit 
a dual description in terms of string theories, defined in a higher-dimensional spacetime. 
In particular, the "holographic dictionary" relating quantities between the two types of 
theories states that, when the number of color charges in the gauge theory tends to infinity, 
the string coupling g s in the dual string theory tends to zero. Under these conditions, loop 
effects on the string side of the correspondence can be neglected, i.e. the theory reduces 
to its "classical string" limit. (In addition, the limit in which the gauge theory is strongly 
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coupled corresponds to the limit in which the string length tends to zero. Hence, if the 
gauge theory has a large number of colors and is strongly coupled, its string dual reduces to 
a classical gravity theory, in an appropriate, curved, higher-dimensional spacetime, opening 
up the possibility of analytical treatment). 

The first explicit version of this conjecture was formulated by studying a system of 
N D3 branes (i.e. 3 + 1-dimensional hyperplanes) in type IIB string theory, defined in a 
ten-dimensional spacetime given by the direct product of a five dimensional anti-de Sitter 
spacetime with a five-dimensional sphere, AdSs xS 5 : it was realized that the low-energy 
dynamics of such a system can be equivalently described in terms of supersymmetric Yang- 
Mills (SYM) theory, with U(iV) gauge group and M = 4 supercharges, and that the two 
theories share the same global symmetries. The parameters of the gauge theory (the number 
of colors N and the 't Hooft coupling A) are mapped to those of the string theory (the string 
coupling g s , and the ratio of the string length l s to the radius of the space R) as: 

X/N = A7rg s , X = (l s /R)~ 4 . (5.1) 

For the J\f = 4 SYM theory, the validity of the gauge/string correspondence is supported 
by many pieces of mathematical evidence, and no counter-examples to it have ever been 
found (although a general proof has not been formulated yet). // the correspondence is 
indeed true, and generic (i.e., if any gauge theory admits a holographic description in 
terms of a string model), then it would be possible to study the strong-coupling regime 
of gauge theories (at least of those with a large number of colors) via simple, "classical" 
calculations on the string theory side. This possibility has motivated a huge theoretical 
effort, aimed at extending the holographic correspondence to theories which are more "QCD- 
like" — including, in particular, theories with matter fields in the fundamental representation 
of the gauge group, exhibiting a non-trivial running coupling, with a spectrum of confined 
states and spontaneous chiral symmetry breaking at low energies. This topic is extensively 
discussed in ref. [14]; in the following, we summarize the main issues relevant for our present 
discussion, and for a comparison with lattice results. 

Before reviewing the holographic approach to meson spectra in the large- N limit, how- 
ever, we would like to warn the reader that, while conceptually extremely interesting, these 
types of computations are not yet at a level of accuracy which enables precise comparison 
with lattice results (or with experimental data). Although the gauge/string correspondence 
has led to dramatic theoretical progress in several aspects of strongly coupled field theories, 
its application to quantitatively address questions of a direct, phenomenological relevance 
is still limited. For these reasons, it should be understood that the theoretical predictions 
discussed in this context have to be taken cum grano salis. 

First of all, it is worth emphasizing that the M = 4 SYM theory, at least at zero or low 
temperatures, is qualitatively very different from QCD: in particular, it is non-confining, 
it does not have a chiral condensate nor a discrete spectrum, it is maximally supersym- 
metric, conformally invariant, and it does not have elementary fields in the fundamental 
representation of the gauge group. One possible way to modify the AdS/CFT correspon- 
dence, in order to obtain a gauge theory with features closer to those of QCD is based on 
the so-called "top-down" approach, i.e., on a deformation of the original AdS/CFT setup 
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with the addition of extra ingredients. In particular, in order to have a gauge theory with 
fields in the fundamental representation of the gauge group, one can modify the dual string 
model by adding a set of nf "flavor" D7 branes [92] (see also, e.g., ref. [93] for a discussion): 
open strings stretching from the original D3 branes to one of these D7 branes correspond 
to fundamental matter fields in the gauge theory, with a mass proportional to the separa- 
tion between the D7 branes and the D3 branes. This reduces the supersymmetry of the 
corresponding gauge theory down to M = 2, and leads to a global U(n/) "flavor" symmetry 
(which is the remnant of the gauge symmetry on the D7-branes). On the other hand, open 
strings that connect two D7 branes are interpreted as the holographic duals of "mesons". In 
this context, "mesons" denote tightly bound states of a quark and its antiquark. 

Assuming that the number rif of D7 branes remains finite, and hence much smaller 
than the number N of D3 branes (which is taken to infinity), one can work in the so- 
called "probe approximation" [92], in which the D7 branes do not affect the geometry of the 
spacetime in which the string theory is defined. It is interesting to note that, in the dual 
gauge theory, this corresponds to a "quenched approximation" (which, historically, has been 
used for a long time also in numerical simulations of lattice QCD [44]), in which quarks 
are propagating in a background generated by the gauge fields only. This approximation 
becomes exact in the 't Hooft limit, in which the 1/N suppression of virtual quark loops 
leads to a unitary, quenched theory, whose dynamics is purely determined by planar gluon 
loops. 

Within the probe approximation in holography, one can easily derive the structure of 
the meson spectrum: for quarks of non-vanishing mass m q (that is, for a non-vanishing 
separation between the D3 and the D7 branes), one finds the spectrum to be [94]: 



MH=4^J^ + 1 f + 2 ) 



(5.2) 



where n denotes the quantum number describing the radial excitation of the meson. Some 
observations are in order: 

1. eq. (5.2) shows that, in the limit of light quarks, the meson mass vanishes — and it 
does so linearly in m q . This is in contrast to the situation in QCD, where, in general, 
in the chiral limit mesons remain massive (due to the existence of a mass gap) or, in 
the case of pseudoscalar mesons which are the Nambu-Goldstone bosons associated 
with the spontaneous breakdown of chiral symmetry, with masses that vanish like 
,fm~ n 



the dependence of M^ n > on n described by eq. (5.2) is different from the one char- 
acteristic of Regge trajectories (M^ a > oc \/n): in principle, the correct functional 
dependence characterizing the spectrum in the large- TV limit could unambiguously be 
identified through lattice calculations. In practice however it turns out that the lat- 
tice determination of even the first few excited states is a technically very demanding 
task (see ref. [95] and references therein for a discussion), hence in this paper we do 
not attempt a systematic computation in this direction, and limit our investigation 
to the first excited state in some channels only; 
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3. at finite quark mass and finite but large 't Hooft coupling A, eq. (5.2) implies that 
the mass of the lightest meson is finite, but parametrically suppressed with respect 
to m g : 

m<°>=47rm tf y|, (5.3) 

and, hence, that the constituent quark and antiquark form a tightly bound state. 

With the caveats mentioned above, it is nevertheless possible to consider a quantitative 
comparison of our results with holographic predictions in the probe approximation. As an 
example, we consider the dependence of the p vector meson mass on the mass of the pion, 
which can be compared with the prediction derived in ref. [93] in a setup involving a 
background with a non-constant dilaton [73]: 



p \ ' -1 + 0.307 
m p (0) 



ITIr, 



m p (0) 



(5.4) 



Our numerical result for the same quantity in the large- N limit, reported in eq. (3.22), is 
remarkably close to the above expectation. 

The inclusion of backreaction effects of the D7 branes on the spacetime geometry makes 
the meson spectrum computation more complicated, but it is still possible to study the 
meson mass dependence for a probe brane in a background which includes curvature effects 
due to the flavor branes. The result reads: 

M 2 oc (n+l)(n + 2)m 2 q hx(m 2 q ) . (5.5) 

Various authors have proposed different constructions, characterized by less super- 
symmetry — see, e.g., refs. [96-100]. In particular, ref. [98] discusses a setup involving a set 
of D4, D8 and D8 branes, with a compactified direction of radius inversely proportional to 
a characteristic Kaluza-Klein scale, and derives a dual description featuring massless pions 
and interesting predictions for other meson masses. According to this construction, the 
ratio between the squared masses of the states with the quantum numbers corresponding 
to the ai(1260) and the p mesons turns out to be: 



2 
m ai(1260) 

771,2 



r~^ 



2.4, (5.6) 



which is compatible with the experimental unquenched N = 3 value 2.5(1) [101]. This 
value can also be compared to our result for this quantity, extrapolated to the chiral and 
large- A limits, which is 3.6(3). 

An even less favorable comparison holds for the ratio of the fundamental and first 
excited states in the vector channel: 

^^4.3. (5.7) 

777~ 
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While the experimentally observed value for this quantity is 3.57(12) [101], our N = oo 
result reads 6.0(4). Regarding the (isovector) scalar channel, the prediction of this model 
for the squared ratio of the lightest state mass to the lightest vector reads: 

2 
m ao(1450) ^ 

to be compared with the real-world value 3 3.61(9) [101]. Our result, extrapolated to van- 
ishing quark mass and for N — > oo, is 8.4(6). 

Finally, we mention that a different type of strategy for holographic studies was pio- 
neered in refs. [102, 103]: it goes under the name of "AdS/QCD", and consists of constructing 
a gravity dual in a curved higher-dimensional (typically: five-dimensional) space, reproduc- 
ing the known features of QCD (a related approach has been discussed in refs. [104-109]). 
Contrary to the constructions mentioned above, here one follows a "bottom-up" approach, 
which would be related to a non-critical string theory setup [110, 111]. An important 
caveat in such constructions, however, is that string corrections may be quantitatively 
non-negligible at finite values of the 't Hooft coupling, hence in this case the gravity ap- 
proximation could be unjustified. However, these effective models appear to capture salient 
features of QCD, including, in particular, confinement and chiral symmetry breaking, and 
typically yield quantitatively rather accurate predictions for certain physical quantities (or 
ratios thereof). The quantitative predictions of the model proposed in ref. [102] were worked 
out using two different methods. One possibility is to fix the values of the free parameters 
of the model, by setting the masses of 7r and p, and the pion decay constant to their physical 
values. This led to the predictions: 

m Ol(1260 ) = 1363 MeV, F p = (329 MeV) 2 (5.9) 

(note that, upon conversion to our notations, this corresponds to f p = 198 MeV). Alter- 
natively, one can perform the best fit for all of the seven parameters simultaneously, which 
results in: 

m^ = 141 MeV, m p = 832 MeV, m ai(1260 ) = 1220 MeV, 

F n = 84 MeV, F p = (353 MeV) 2 (5.10) 

(which, in our conventions, would correspond to f p = 212 MeV). 

In principle, these values can also be compared to our results extrapolated to the 
large- A limit. Note, however, that, given that the parameters of this model involve input 
from experimental data, the values obtained for N — > oo are not necessarily expected 
to be in better agreement than those for N = 3. However, our lattice computations at 
different values of iV can provide helpful insights into the consistency of the model: since 



3 Strictly speaking, the lightest experimentally measured state with quantum numbers 7 (J ) = 
1~(0 ++ ) is the ao(980), rather than the ao(1450). However, the correct identification of scalar mesons 
is a particularly challenging problem, and the identification of the ao(980) as a genuine meson is somewhat 
controversial: most likely, the ao(980) wave function has a large KK component [101], and hence this state 
could be interpreted as a two-meson resonance or a tetraquark. 
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the holographic construction is a gravitational one, i.e. it is based on the approximation 
of an infinite number of colors, the lattice results can reveal the quantitative impact of 
corrections due to the finiteness of N, and therefore provide a non-trivial test of the validity 
of the model. Our results reveal that, in most cases, the finite-iV corrections evaluated at 
N = 3 amount for relative corrections well below 10%. 

5.2 Chiral perturbation theory 

As is well-known, chiral perturbation theory [112-117] is an effective low-energy theory de- 
scribing the dynamics of the lightest mesons in QCD. It relies on the parametric separation 
of the chiral symmetry breaking scale ~ 4-7ri ? 7I - and the (nearly) zero mass of the (pseudo-) 
Goldstone bosons. In the case of a QCD-like theory with nf flavors of light quarks, xPT 
describes the fields associated with the light mesons, in terms of the components of a field 
taking values in U(n/). Its dynamics is governed by an effective Lagrangian (constrained by 
the symmetries of the theory), which can be organized in a systematic expansion, according 
to the number of derivatives and of factors involving a possible explicit mass term, in which 
the coefficients of the different terms are low-energy constants (LECs), whose numerical val- 
ues can be fixed using phenomenological input, and compared with the expectations from 
large- N counting rules. In particular, inspection of the terms contributing to the lowest 
order shows that the effective Lagrangian is proportional to the square of the pion decay 
constant, i.e. to N. Essentially, this implies that, for N —> oo, the effective theory for light 
mesons becomes exact at tree level. However, we should point out that the construction of 
%PT for large- N QCD (and the constraints that can be derived from it) is a topic which 
involves some subtleties: the interested reader can find a clear exposition of this subject in 
ref. [118]. In addition, there exist studies of large- iV xPT f° r t ne baryonic sector [119], too. 
While a systematic comparison of large- N xPT predictions with lattice results (in- 
cluding, in particular, the study of the iV-scaling of various LECs) is a task that would 
go way beyond the scope of the present work, for our present purposes it is worthwhile 
mentioning refs. [84-86], in which the iV-dependence of the masses for some of the mesons 
presented here was investigated in xPT. These studies address the full theory with sea 
quarks. In this case meson masses tended to increase with the number of colors relative to 
F = W3/NFn(0), as opposed to our quenched results. The present work and xPT results 
can be compared by studying ratios of quantities defined at N = oo, where both theories 
are quenched. Comparing our results shown in table 4 with those in refs. [84-86] we find a 
discrepancy of order 30%, however the systematics of their approach are hard to estimate, 
in particular because experimental data are only available for N = 3. In the future our 
N = oo results could be used to constrain effective field theory predictions. 

6 Conclusions 

We have computed decay constants as well as the ground and first excited state masses of 
mesons in the large- N limit of QCD by lattice simulation of the N = 2, 3, 4, 5, 6, 7 and 17 
quenched theories. In all channels but the scalar the dependence on the number of colors 
for iV > 3 is mild. The corrections are well parameterised by an expansion in 1/N 2 . 
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Extrapolating our quenched results to the large- A limit (which is unitary and the same 
as that of QCD with sea quarks), we determine the meson spectrum in the chiral limit as 
well as at physical light and strange quark masses. We find the scalar to be about 1.5 times 
heavier than the vector particle at quark masses smaller than that of the strange quark. 
This is of particular relevance to the phenomenology of scalar mesons [84-86] . 

We also clarified a discrepancy among previous studies: our results for the p meson 
mass for SU(17) are compatible with the large- iV extrapolation of studies carried out for 
smaller values of N [45-48], and are in contrast to the findings obtained in ref. [49] at 
N = 17, 19 using different techniques. This suggests that the disagreement with the rest of 
the literature may be due to some technical aspect of ref. [49] . 

Finally, we compared our numerical results with analytical predictions, including, in 
particular, some of those derived from holographic models, finding qualitative (and semi- 
quantitative) agreement. Our results for the masses of various states (expressed in units 
of the p meson mass) exhibit a systematic tendency towards values which are larger than 
those obtained from holographic computations. Expressed in units of the string tension (or 
of the pion decay constant normalized by vN), we find meson masses to decrease with an 
increasing number of colors. This may be different in QCD with sea quarks. 
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A Renormalization constants 

The perturbative expansions for the non-singlet renormalization constants Za and Zy are 
known to order A 2 . As discussed above, in this work we considered two definitions of the 
quark mass, namely: 

1/1 1 \ 1 ,. ^ 

m c , (A.l) 



2a \k k c J 2an 
and the definition through the axial Ward identity: 

_ (0|VUx)|7t) 

mPCAC - 2<0|P(*)|7T> ' (A ' 2) 

where A^ = q^^^q and P = q'j^q are the non-singlet axial and pseudoscalar local currents 
for a quark of mass m q . 
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Converting the lattice results to the MS scheme at a scale [i = n/a amounts to: 

m Ms( 7T / a ) = ~^~ m i = ^- m PCAc, (A. 3) 



with the renormalization constants [77, 78]: 



Zs = l- 0.515360027(4) ( 1 - — ^ J — - ( 0.6155(2 



1 A A f n „ KKfn , 0.7325(2) \ / l\(\\ 



2 






N 2 J V N 2 J \4ttJ ' 

(A.4) 



* - . - o_ (4) (x - i,) A _ («, - S-Jffl) (, - i) (A) , 

(A.5) 
Z v = l- 0.8203561429(3) M - ^ ) ^ " ( 1.01790(5) - °^^^1 j f . \ ( A 



A 2 / V4vr 

(A.6) 



Za = 1 - 0.4693595879(2) ( 1 - — ^ j ( 0.09173(5) - N 2 ) f 1 



1 A A fn nM »of K \ 0.15084(2)^ (^ X \( X 

N 2 ) V47T 



2 



(A.7) 

where A = 2N 2 //3 is the 't Hooft coupling in the lattice scheme. Z$ and Zp above are 
defined at the lattice cut-off scale [i = n/a, while Zy and Za are scale-independent. The 
scale dependence also cancels from the ratio Zs/Zp. For the sake of completeness, we also 
quote the expansion of m c = l/(2aK c ) [120, 121]: 

1 \ A 

2 



m r = 4 - 2.046522156925001893 1 , , 

N 2 J 4tt 



^mm-""-®)^)®. (A ,, 



As is well-known, lattice perturbation theory converges slowly, but the convergence can 
usually be boosted, using an improved definition of the expansion parameter, e.g., Aj^g(7r/a), 
taking the latter value from a lattice computation. Alternatively, one can replace the bare 
lattice coupling with a value computed from a lattice observable [122, 123]. 

Expanding the average plaquette 

U D = 1 - ciA - c 2 A 2 + . . . (A.9) 



leads to one- and two-loop improved couplings of the form 



where [124, 125]: 





»S'= 


-* 


lnU n , 






4 2> = 


*S> 


- did 2 (In Up) -, 




>h 


1 

Cl 


-( 


}-^)' 




>h 


1 

= 2 + 


C2 

cf 


0. 826844 A 2 - 


1 


N 2 -l 





(A.10) 
(A.ll) 

(A.12) 
(A.13) 
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N 


P 


U a 


A 


A (1) 
A E 


Ag> 


2 


2.4645 


0.644718(3) 


3.246 


4.682 


3.101 


3 


6.0175 


0.595595(3) 


2.991 


4.664 


2.718 


4 


11.0280 


0.578794(2) 


2.902 


4.666 


2.586 


5 


17.5350 


0.573069(2) 


2.851 


4.640 


2.522 


6 


25.4520 


0.568682(2) 


2.829 


4.644 


2.490 


7 


34.8343 


0.566504(2) 


2.813 


4.641 


2.471 


17 


208.45 


0.562729(4) 


2.773 


4.616 


2.423 



Table 5. Plaquette values, and 't Hooft couplings: one- and two- loop improved (E). 



The resulting values are displayed in table 5. To first and second order accuracy in pertur- 
bation theory A^, and X E agree with A, respectively. 

In order to obtain improved 0(A) results, we replace A by A^. in eqs. (A.4)-(A.7) 

above, while for improved 0(A 2 ) renormalization constants we replace A by A^. and A 2 by 

(l - ) 2 
X E , for consistency to this order in perturbation theory. 

In table 6, we display the resulting perturbative renormalization constants and critical 
hopping parameters. Our non-perturbative results on k c can be found in table 8. The one- 
loop improved Z4- and Zy-factors seem to be closest to the corresponding non-perturbative 
results [68, 76] (see tables 6,7) obtained for N = 3 at a similar coupling (A = 3.0 vs. our 
A = 2.991). In the absence of non-perturbative results for A 7^ 3, we use these "1-loop E" 
Za- and Zy-factors also for the other SU(iV) groups, opting for a conservative 8% estimate 
of their relative uncertainties. Note that the differences at N = 3 to the non-perturbative 
results are 7.5 % and 6.5 %, respectively. 

The quantity Zp/(ZaZs), needed to relate the quark to the PCAC mass, is listed as 
well in table 6. We find a huge difference between our perturbative expectations and the 
non-perturbative results for N = 3 of refs. [68, 76], obtained at A = 3.0 and listed in table 7. 
Therefore, we fix the latter ratio non-perturbatively from the fit eq. (3.7). 



32 





N 


1-loop 


2-loop 


1-loop E 


2-loop E 


Zs 


2 


0.9 


0.879 


0.856 


0.86 


3 


0.891 


0.864 


0.83 


0.836 




4 


0.888 


0.86 


0.821 


0.827 




5 


0.888 


0.859 


0.817 


0.824 




6 


0.887 


0.858 


0.815 


0.822 




7 


0.887 


0.857 


0.814 


0.821 




17 


0.887 


0.857 


0.811 


0.819 


Z P 


2 


0.887 


0.857 


0.838 


0.829 


3 


0.877 


0.84 


0.808 


0.798 




4 


0.874 


0.835 


0.798 


0.787 




5 


0.873 


0.834 


0.794 


0.783 




6 


0.873 


0.833 


0.791 


0.78 




7 


0.873 


0.832 


0.79 


0.778 




17 


0.872 


0.832 


0.787 


0.776 


Zv 


2 


0.841 


0.801 


0.771 


0.764 


3 


0.826 


0.78 


0.729 


0.729 




4 


0.822 


0.774 


0.714 


0.717 




5 


0.821 


0.773 


0.709 


0.713 




6 


0.82 


0.771 


0.705 


0.71 




7 


0.82 


0.771 


0.703 


0.708 




17 


0.82 


0.77 


0.7 


0.706 


Z A 


2 


0.909 


0.906 


0.869 


0.907 


3 


0.901 


0.897 


0.845 


0.901 




4 


0.898 


0.894 


0.837 


0.899 




5 


0.898 


0.894 


0.834 


0.898 




6 


0.897 


0.893 


0.831 


0.898 




7 


0.897 


0.893 


0.83 


0.898 




17 


0.897 


0.892 


0.828 


0.898 




2 


0.13875 


0.14295 


0.14585 


0.14702 


K c 


3 


0.14018 


0.14526 


0.15038 


0.15137 




4 


0.14057 


0.145912 


0.15209 


0.15298 




5 


0.14068 


0.14609 


0.15269 


0.15359 




6 


0.14076 


0.14623 


0.15316 


0.15396 




7 


0.14080 


0.14629 


0.15339 


0.15417 




17 


0.14085 


0.14637 


0.15380 


0.15453 


Z P 


2 


1.085 


1.076 


1.126 


1.063 


Za-Zs 


3 


1.093 


1.084 


1.152 


1.06 




4 


1.095 


1.086 


1.162 


1.058 




5 


1.096 


1.087 


1.166 


1.057 




6 


1.097 


1.087 


1.168 


1.057 




7 


1.097 


1.087 


1.17 


1.056 




17 


1.097 


1.087 


1.172 


1.056 



Table 6. Perturbative Z-factors and n c : one-loop, two-loop, bare and improved (E). 
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Ref. 

w 



Zs 
0.68(1) 



Z P 
0.45(6) 



Z A 
0.81(1) 



Zy 
0.71(2) 



Zp/(ZaZs) 
0.82(11) 



[76] 



0.7718(16) 



0.4934(24) 0.7821(9) 



0.6833(8) 0.8174(44) 



Table 7. Non-perturbative Z-factors from refs. [68, 76] obtained for N = 3 at A = 3.0 (close to 
our A = 2.991). 



B Additional tables and figures 

This appendix includes additional tables and figures. 
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Figure 15. Fit of the a,\ mass to: m ai /^fcr = A + B ■ m PCAC /y/<T (left) and 1/N 2 fit of the 
parameters A and B (right). 
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Figure 16. Same as figure 15, but for the mass of the b\ state. 
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Figure 17. Same as figure 15, but for the mass of the tt* state. 
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Figure 18. Same as figure 15, but for the mass of the p* state. 
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Figure 19. Same as figure 15, but for the mass of the a* state. 
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Figure 20. Same as figure 15, but for the mass of the b\ state. 




Figure 21. Same as figure 15, but for the mass of the q,q state. 
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N 


Zp/(ZaZs) 


A 


K c 


2 


0.6831(75) 


0.972(74) 


0.152880(22)(81) 


3 


0.7519(42) 


0.739(38) 


0.156670(8)(86) 


4 


0.7812(36) 


0.607(32) 


0.158218(7)(87) 


5 


0.8050(25) 


0.483(21) 


0.158596(5) (87) 


6 


0.8098(22) 


0.462(18) 


0.159103(4) (87) 


7 


0.8133(19) 


0.447(16) 


0.159326(4)(87) 


17 


0.8286(55) 


0.370(42) 


0.159590(12)(87) 



Table 8. Fit parameters of the PCAC mass. 



N 


A Kc /A/3 


A/3/A(a^) 


5(ay/a) 


5k c 


2 


0.0588172 


1.37097 


0.001 


0.0000806 


3 


0.033977 


2.54391 


0.001 


0.0000864 


4 


0.019335 


4.4773 


0.001 


0.0000866 


6 


0.00891258 


9.78578 


0.001 


0.0000872 



Table 9. Systematic errors of k c and its dependencies, according to eq. (3.6). 



K, 


Nf xN t 


mpcAc/y^a 


"v/v^ 


rrip/y/a 


maJy/fJ 


mbJVcf 


mao/y/cf 


0.14581 


16 3 x 32 

24 3 x 48 
32 3 x 64 


0.5786(11) 
0.57866(90) 
0.5776(12) 


2.7042(85) 
2.6926(56) 
2.6835(58) 


2.931(13) 
2.9110(82) 
2.9122(93) 


4.038(50) 
4.136(64) 
4.133(76) 


4.043(49) 
4.115(58) 
4.125(70) 


3.952(59) 
4.110(52) 
4.085(76) 


0.14827 


16 3 x 32 

24 3 x 48 
32 3 x 64 


0.3551(12) 
0.35564(84) 
0.3548(11) 


2.131(10) 
2.1116(61) 
2.0968(59) 


2.464(18) 
2.438(10) 
2.431(11) 


3.622(72) 
3.730(79) 
3.707(92) 


3.602(69) 
3.743(75) 
3.679(78) 


3.565(75) 
3.745(72) 
3.71(11) 


0.15008 


16 3 x 32 

24 3 x 48 
32 3 x 64 


0.2054(14) 
0.20732(83) 
0.2065(10) 


1.664(12) 
1.6322(72) 
1.6138(63) 


2.134(24) 
2.104(16) 
2.086(18) 


3.33(11) 
3.43(11) 
3.40(11) 


3.29(10) 
3.529(92) 
3.353(94) 


3.49(11) 

3.65(12) 
3.57(20) 


0.15096 


16 3 x 32 

24 3 x 48 
32 3 x 64 


0.1363(17) 
0.14047(88) 
0.13964(99) 


1.400(15) 
1.3617(89) 
1.3444(68) 


1.977(29) 
1.947(23) 
1.933(26) 


3.18(16) 
3.27(13) 
3.25(14) 


3.14(15) 
3.50(12) 
3.22(12) 


3.97(22) 
3.96(15) 
3.68(39) 


0.151959 


24 3 x 48 
32 3 x 64 


0.0661(15) 
0.0660(11) 


0.960(15) 
0.9558(100) 


1.818(32) 
1.813(40) 


3.24(14) 
3.17(16) 


3.75(20) 
3.38(20) 


1.97(36) 
3.67(16) 


0.152496 


24 3 x 48 
32 3 x 64 


0.0236(23) 
0.0186(18) 


0.639(39) 
0.538(30) 


1.812(71) 
1.730(60) 


3.29(39) 
3.18(15) 


5.1(1.3) 
3.64(26) 


2.88(32) 



Table 10. Ground state meson masses of SU(2) gauge theory. 
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K 


iV 3 x N t 


"T-PCAc/^ 


in n / 'y/a 


mp/y/a 


m ai /y/a 


nibjy/a 


m ao /\/cr 


0.15002 


16 3 x 32 
24 3 x 48 
32 3 x 64 


0.5514(11) 
0.55154(55) 
0.55138(58) 


2.6597(73) 
2.6528(31) 
2.6513(33) 


2.915(11) 
2.9066(47) 
2.9033(58) 


4.081(35) 
4.068(36) 
4.117(57) 


4.090(33) 
4.071(44) 
4.129(54) 


3.938(37) 
3.930(27) 
3.982(43) 


0.1522 


16 3 x 32 
24 3 x 48 
32 3 x 64 


0.3555(11) 
0.35558(52) 
0.35558(51) 


2.1186(86) 
2.1109(32) 
2.1111(32) 


2.476(13) 
2.4642(57) 
2.4631(70) 


3.678(42) 
3.684(43) 
3.739(68) 


3.684(42) 
3.674(42) 
3.769(65) 


3.512(50) 
3.518(37) 
3.577(67) 


0.1538 


16 3 x 32 
24 3 x 48 
32 3 x 64 


0.2209(12) 
0.22105(49) 
0.22132(47) 


1.669(10) 
1.6606(36) 
1.6626(34) 


2.156(16) 
2.1362(76) 
2.1399(92) 


3.390(55) 
3.400(55) 
3.464(86) 


3.415(56) 
3.386(53) 
3.529(83) 


3.266(83) 
3.248(59) 
3.299(91) 


0.15458 


16 3 x 32 
24 3 x 48 
32 3 x 64 


0.1582(13) 
0.15864(48) 
0.15907(44) 


1.418(12) 
1.4099(40) 
1.4130(35) 


2.006(20) 
1.9781(98) 
1.979(13) 


3.260(68) 
3.259(64) 
3.32(10) 


3.323(73) 

3.242(64) 
3.422(98) 


3.25(12) 
3.169(90) 
3.163(98) 


0.155638 


24 3 x 48 
32 3 x 64 


0.07741(50) 
0.07818(42) 


0.9993(56) 
1.0030(39) 


1.788(19) 
1.796(24) 


3.134(62) 
3.064(87) 


3.148(80) 
3.26(10) 


3.51(21) 
3.20(26) 


0.15613 


24 3 x 48 
32 3 x 64 


0.03958(69) 
0.04128(43) 


0.7306(86) 
0.7411(54) 


1.710(25) 
1.749(26) 


3.032(93) 
2.88(11) 


3.15(18) 
3.23(17) 


2.51(23) 
3.07(61) 



Table 11. SU(3) ground state meson masses. 



K 


iV 3 x N t 


mp CAC /y/a 


m^/y/a 


mp/y/a 


majy/cr 


m bl /y/a 


mao/y/a 


0.15184 


16 3 x 32 

24 3 x 48 


0.52889(72) 
0.52863(38) 


2.6147(71) 
2.6168(24) 


2.882(11) 
2.8921(38) 


3.999(41) 
4.070(28) 


4.050(44) 
4.104(24) 


3.855(39) 
3.914(22) 


0.154 


16 3 x 32 

24 3 x 48 


0.33716(70) 
0.33710(36) 


2.0641(83) 
2.0665(24) 


2.437(16) 
2.4513(46) 


3.575(56) 
3.656(29) 


3.639(60) 
3.691(30) 


3.406(56) 
3.482(29) 


0.15559 


16 3 x 32 

24 3 x 48 


0.20413(73) 
0.20423(35) 


1.5965(97) 
1.5993(25) 


2.106(25) 
2.1227(64) 


3.235(80) 
3.343(37) 


3.331(91) 
3.363(40) 


3.137(65) 
3.184(47) 


0.15635 


16 3 x 32 

24 3 x 48 


0.14303(77) 
0.14332(35) 


1.337(11) 
1.3390(26) 


1.954(31) 
1.9651(85) 


3.04(10) 
3.198(46) 


3.19(13) 
3.198(49) 


3.03(11) 
3.095(73) 


0.157173 


24 3 x 48 


0.07935(36) 


1.0001(32) 


1.802(13) 


3.088(80) 


3.009(87) 


3.26(12) 


0.15764 


24 3 x 48 


0.04350(39) 


0.7498(40) 


1.709(20) 


3.03(12) 


3.00(16) 


3.34(26) 



Table 12. SU(4) ground state meson masses. 
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K 


iV 3 x N t 


m PCACI \f& 


fn-wl ' \[o 


ITlp/y/cf 


maJVa 


m^/y/a 


rriajx/a 


0.15205 


16 3 x 32 

24 3 x 48 


0.55192(54) 
0.55112(33) 


2.6736(37) 
2.6680(18) 


2.9321(55) 
2.9331(28) 


4.070(28) 
4.080(20) 


4.114(31) 
4.127(23) 


3.872(30) 
3.904(20) 


0.15426 


16 3 x 32 

24 3 x 48 


0.35403(54) 
0.35314(31) 


2.1156(43) 
2.1065(20) 


2.4774(72) 
2.4769(36) 


3.671(36) 

3.658(25) 


3.711(43) 
3.708(25) 


3.420(47) 
3.446(27) 


0.15592 


16 3 x 32 

24 3 x 48 


0.21349(56) 
0.21247(29) 


1.6322(54) 
1.6192(22) 


2.1292(100) 
2.1263(54) 


3.374(54) 
3.350(33) 


3.402(67) 
3.411(32) 


3.095(61) 
3.115(46) 


0.15658 


16 3 x 32 

24 3 x 48 


0.15960(58) 
0.15852(29) 


1.4102(61) 
1.3954(24) 


1.989(12) 
1.9851(69) 


3.257(69) 
3.238(43) 


3.276(90) 
3.304(37) 


3.021(98) 
3.006(66) 


0.157548 


24 3 x 48 


0.08158(29) 


1.0017(29) 


1.779(13) 


3.113(60) 


3.201(72) 


2.96(12) 


0.158355 


24 3 x 48 


0.01833(38) 


0.4869(68) 


1.655(32) 


3.076(64) 


3.49(14) 


2.80(35) 



Table 13. SU(5) ground state meson masses. 



K 


TV 3 x N t 


r7l PCAc/\/ Cr 


fn-wl ' \[o 


rrip/y/a 


maJVa 


m bl /y/a 


mao/Vcr 


0.15264 


16 3 x 32 

24 3 x 48 


0.54200(48) 
0.54185(26) 


2.6579(35) 
2.6574(16) 


2.9299(51) 
2.9279(25) 


4.064(29) 
4.069(44) 


4.091(27) 
4.104(17) 


3.910(29) 
3.863(17) 


0.15479 


16 3 x 32 

24 3 x 48 


0.35056(49) 
0.35042(24) 


2.1079(42) 
2.1086(17) 


2.4868(67) 
2.4832(30) 


3.630(40) 
3.654(21) 


3.664(38) 
3.694(20) 


3.456(42) 
3.388(22) 


0.15636 


16 3 x 32 

24 3 x 48 


0.21769(51) 
0.21774(23) 


1.6436(50) 
1.6464(18) 


2.1549(90) 
2.1504(42) 


3.355(39) 
3.338(24) 


3.392(35) 
3.387(24) 


3.127(72) 
3.018(31) 


0.15712 


16 3 x 32 

24 3 x 48 


0.15545(53) 
0.15570(22) 


1.3820(57) 
1.3867(20) 


1.990(12) 
1.9864(54) 


3.206(51) 
3.180(27) 


3.232(47) 
3.235(29) 


3.00(11) 

2.834(42) 


0.158051 


24 3 x 48 


0.08177(22) 


1.0013(25) 


1.7819(100) 


2.944(60) 


3.001(76) 


2.687(65) 


0.158845 


24 3 x 48 


0.01948(50) 


0.4989(55) 


1.569(36) 


2.871(65) 


3.46(16) 


2.39(25) 



Table 14. SU(6) ground state meson masses. 



K 


TV 3 x N t 


mpcAc/y/v 


m^/y/a 


mp/y/a 


rriajy/a 


m bl /y/a 


mao/y/a 


0.152816 


16 3 x 32 

24 3 x 48 


0.54657(39) 
0.54604(22) 


2.6767(31) 
2.6727(13) 


2.9489(47) 
2.9461(20) 


4.105(26) 
4.117(13) 


4.117(27) 
4.140(14) 


3.913(36) 
3.913(15) 


0.154967 


16 3 x 32 

24 3 x 48 


0.35511(39) 
0.35455(21) 


2.1322(34) 
2.1242(14) 


2.5100(55) 
2.5009(25) 


3.689(32) 
3.709(16) 


3.712(33) 
3.730(17) 


3.464(25) 
3.441(19) 


0.156547 


16 3 x 32 

24 3 x 48 


0.22147(41) 
0.22093(20) 


1.6689(40) 
1.6599(16) 


2.1781(75) 
2.1631(36) 


3.356(46) 
3.405(20) 


3.412(51) 

3.427(22) 


3.147(30) 
3.070(28) 


0.157339 


16 3 x 32 

24 3 x 48 


0.15669(42) 
0.15619(20) 


1.3983(45) 
1.3892(17) 


2.0096(98) 
1.9895(48) 


3.171(61) 
3.251(26) 


3.262(72) 
3.272(27) 


3.013(47) 

2.874(38) 


0.158273 


24 3 x 48 


0.08183(20) 


1.0007(19) 


1.7785(81) 


3.076(54) 


3.026(65) 


2.704(62) 


0.159062 


24 3 x 48 


0.02034(24) 


0.4960(40) 


1.594(28) 


2.917(54) 


2.930(73) 


3.09(30) 



Table 15. SU(7) ground state meson masses. 
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K 


iV 3 x N t 


W-PCAc/xA 7 


fTl-xl \fo 


rrip/^/a 


m ai /y/a 


m^lsfa 


m ao /yfa 


0.15298 


12 3 x 24 


0.55913(42) 


2.7030(56) 


2.9633(77) 


4.122(23) 


4.170(24) 


3.917(24) 


0.15521 


12 3 x 24 


0.35983(44) 


2.1356(69) 


2.499(10) 


3.681(33) 


3.733(34) 


3.423(37) 


0.15684 


12 3 x 24 


0.22114(46) 


1.6545(86) 


2.148(14) 


3.344(46) 


3.410(50) 


3.016(60) 


0.15755 


12 3 x 24 


0.16257(47) 


1.4104(99) 


1.992(17) 


3.191(54) 


3.277(63) 


2.812(84) 


0.158531 


12 3 x 24 


0.08390(68) 


0.998(12) 


1.776(22) 


3.019(70) 


3.161(97) 


2.59(16) 


0.15931 


12 3 x 24 


0.0212(11) 


0.474(38) 


1.580(66) 


3.23(34) 


3.20(29) 


3.3(11) 



Table 16. SU(17) ground state meson masses. 



K 


N 3 x N t 


m-x* J yfo 


m p */^/a 


fn a \l\fo 


m h */^fa 


™> a *ly/a 


0.14581 


16 3 x 32 

24 3 x 48 
32 3 x 64 


4.56(11) 
4.74(15) 
4.69(13) 


4.67(13) 

4.586(77) 
4.86(14) 


5.56(12) 
5.59(19) 
5.62(12) 


5.61(11) 
5.73(20) 
5.84(14) 


5.51(13) 
5.36(12) 
5.34(14) 


0.14827 


16 3 x 32 

24 3 x 48 
32 3 x 64 


4.16(15) 
4.27(19) 
4.28(16) 


4.35(16) 
4.229(90) 
4.52(17) 


5.24(18) 
5.11(32) 
5.26(15) 


5.32(14) 

5.47(28) 
5.61(18) 


5.28(21) 

4.87(18) 
4.87(18) 


0.15008 


16 3 x 32 

24 3 x 48 
32 3 x 64 


3.85(20) 
3.96(27) 
3.95(19) 


4.20(19) 

4.02(10) 
4.27(18) 


4.96(27) 
4.82(23) 
5.06(21) 


5.17(19) 
5.26(30) 
5.31(22) 


4.92(36) 
4.54(34) 
4.59(25) 


0.15096 


16 3 x 32 

24 3 x 48 
32 3 x 64 


3.71(27) 
3.86(30) 
3.75(23) 


4.16(21) 

3.95(11) 
4.14(20) 


4.73(31) 
4.68(29) 
5.00(27) 


5.08(24) 
5.5(1.4) 
5.50(34) 


5.09(67) 

4.87(64) 
4.51(35) 



Table 17. First excited state meson masses of SU(2) gauge theory. 



K 


NfxN t 


m^/yfci 


m p +/yfo 


m a */y/a 


m b */y/a 


m a *Jy/a 


0.15002 


16 3 x 32 

24 3 x 48 
32 3 x 64 


4.599(78) 
4.435(86) 
4.417(82) 


4.73(13) 
4.549(95) 
4.504(79) 


5.666(80) 
5.617(85) 
5.39(11) 


5.821(82) 
5.639(97) 
5.51(13) 


5.535(95) 
5.519(62) 
5.134(96) 


0.1522 


16 3 x 32 
24 3 x 48 
32 3 x 64 


4.239(98) 
4.05(11) 
4.07(11) 


4.39(16) 

4.21(12) 

4.195(100) 


5.43(11) 

5.42(12) 
5.11(16) 


5.66(11) 
5.35(13) 
5.21(16) 


5.26(14) 
5.270(92) 
4.70(14) 


0.1538 


16 3 x 32 
24 3 x 48 
32 3 x 64 


3.99(13) 
3.77(14) 
3.81(17) 


4.18(18) 
3.99(14) 
4.01(12) 


5.33(16) 
5.26(11) 
4.96(16) 


5.70(14) 
5.20(19) 
5.02(21) 


5.12(24) 
5.13(14) 
4.34(20) 


0.15458 


16 3 x 32 
24 3 x 48 
32 3 x 64 


3.88(15) 
3.65(18) 
3.68(23) 


4.11(20) 
3.89(15) 
3.95(14) 


5.35(20) 
5.23(15) 
4.86(20) 


5.82(15) 
5.16(14) 
4.96(26) 


5.22(37) 
5.09(20) 
4.44(18) 



Table 18. SU(3) first excited state meson masses. 
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K 


N% x N t 


m^/y/a 


m p */y/a 


m a */y/a 


m h */y/a 


m a*J\fc 


0.15184 


16 3 x 32 
24 3 x 48 


4.84(11) 
4.576(58) 


4.67(17) 
4.705(65) 


5.65(17) 
5.507(80) 


5.72(21) 
5.531(82) 


5.51(17) 
5.418(48) 


0.154 


16 3 x 32 

24 3 x 48 


4.57(13) 
4.263(74) 


4.39(21) 

4.427(82) 


5.35(22) 
5.19(11) 


5.45(27) 
5.14(11) 


5.20(22) 
5.082(66) 


0.15559 


16 3 x 32 

24 3 x 48 


4.37(16) 
4.07(12) 


4.23(26) 
4.26(10) 


5.08(27) 
4.96(16) 


5.16(31) 

4.82(14) 


4.92(32) 
4.851(94) 


0.15635 


16 3 x 32 

24 3 x 48 


4.26(20) 
4.02(26) 


4.17(30) 
4.20(12) 


4.93(32) 
4.85(20) 


5.00(36) 
4.64(18) 


4.71(43) 
4.75(12) 



Table 19. SU(4) first excited state meson masses. 



K 


iV 3 x N t 


m -K* 1 \[0 


m p */y/a 


mai/y/o 


m b */y/a 


m a */y/a 


0.15205 


16 3 x 32 
24 3 x 48 


4.646(62) 
4.505(47) 


4.795(71) 
4.646(52) 


5.586(89) 
5.471(58) 


5.671(88) 
5.624(67) 


5.36(10) 
5.344(37) 


0.15426 


16 3 x 32 

24 3 x 48 


4.282(79) 
4.146(63) 


4.484(88) 
4.336(68) 


5.32(17) 

5.155(77) 


5.38(13) 
5.300(90) 


4.91(17) 
4.981(48) 


0.15592 


16 3 x 32 

24 3 x 48 


4.00(10) 
3.882(90) 


4.28(10) 
4.125(87) 


5.04(12) 
4.96(11) 


5.18(12) 
5.08(12) 


4.41(25) 
4.754(66) 


0.15658 


16 3 x 32 

24 3 x 48 


3.88(12) 
3.79(13) 


4.20(11) 
4.05(10) 


4.93(15) 
4.89(13) 


5.09(15) 
5.00(15) 


4.10(30) 
4.696(82) 



Table 20. SU(5) first excited state meson masses. 



K 


N% x N t 


m^/y/a 


m p */y/a 


m a * jyfo 


m b */yfa 


m^/V? 


0.15264 


16 3 x 32 
24 3 x 48 


4.494(62) 
4.503(40) 


4.644(72) 
4.634(45) 


5.519(52) 

5.557(48) 


5.674(58) 
5.623(54) 


5.22(16) 
5.396(31) 


0.15479 


16 3 x 32 

24 3 x 48 


4.131(78) 
4.136(53) 


4.329(91) 
4.298(64) 


5.130(73) 
5.231(68) 


5.291(80) 
5.264(74) 


4.83(24) 
5.028(44) 


0.15636 


16 3 x 32 

24 3 x 48 


3.90(10) 
3.874(70) 


4.14(12) 
4.059(92) 


4.846(99) 
4.987(98) 


5.00(11) 
5.00(10) 


4.59(12) 
4.759(65) 


0.15712 


16 3 x 32 

24 3 x 48 


3.82(13) 
3.739(75) 


4.08(14) 
3.94(11) 


4.71(12) 
4.85(13) 


4.85(13) 
4.86(13) 


4.17(28) 
4.633(85) 



Table 21. SU(6) first excited state meson masses. 
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K 


iV 3 x N t 


"V* / \fv 


m p *lsfd 


Tnaxl\[o 


m b */y/a 


m a * Q /Vo- 


0.152816 


16 3 x 32 

24 3 x 48 


4.588(38) 
4.520(36) 


4.743(40) 
4.666(39) 


5.630(43) 
5.585(46) 


5.701(52) 
5.694(59) 


5.386(50) 

5.435(27) 


0.154967 


16 3 x 32 

24 3 x 48 


4.181(49) 
4.118(44) 


4.356(53) 
4.315(49) 


5.297(57) 
5.221(69) 


5.431(57) 
5.338(69) 


5.024(63) 
5.075(36) 


0.156547 


16 3 x 32 
24 3 x 48 


3.928(62) 
3.803(55) 


4.143(64) 
4.061(57) 


5.075(79) 
4.947(83) 


5.190(81) 
5.085(94) 


4.753(92) 
4.839(54) 


0.157339 


16 3 x 32 

24 3 x 48 


3.814(75) 
3.630(67) 


4.054(72) 
3.937(63) 


4.972(100) 
4.82(10) 


5.07(10) 
4.98(12) 


4.63(12) 
4.761(89) 



Table 22. SU(7) first excited state meson masses. 



K 


Nf x N t 


m^/y/a 


m p */y/a 


m a */y/a 


m h */y/a 


ma+JVv 


0.15298 


12 3 x 24 


4.530(48) 


4.660(59) 


5.56(10) 


5.625(85) 


5.344(53) 


0.15521 


12 3 x 24 


4.137(62) 


4.292(76) 


5.199(64) 


5.29(14) 


4.949(81) 


0.15684 


12 3 x 24 


3.865(82) 


4.033(88) 


4.976(86) 


5.07(20) 


4.67(12) 


0.15755 


12 3 x 24 


3.76(10) 


3.928(93) 


4.90(10) 


4.99(25) 


4.57(18) 



Table 23. SU(17) first excited state meson masses. 



State X 


A x ,i 


A X ,2 


Bx,i 


Bx,2 


a 


2.402(34) 


4.25(62) 


2.721(53) 


-6.84(96) 


111 


2.860(21) 


0.84(36) 


2.289(35) 


-2.02(61) 


6i 


2.901(23) 


1.07(40) 


2.273(38) 


-2.83(72) 


TV* 


3.392(57) 


1.0(1.1) 


2.044(80) 


-1.2(1.6) 


P* 


3.696(54) 


0.23(55) 


1.782(67) 


-1.30(54) 


a* 


4.356(65) 


1.8(1.4) 


1.902(98) 


-2.9(2.1) 


a* 


4.587(75) 


1.2(1.2) 


1.76(12) 


-2.1(19) 


b\ 


4.609(99) 


1.7(1.5) 


1.87(15) 


-2.5(2.2) 



Table 24. Coefficients for the expansion of the particle masses as mx/y/v = (Ay,i + A-X.2/N 2 ) + 
(B x ,i + B xa /N 2 ) mpcAc/V^- 



K 


N^ x N t 


F^/^ 


F n /V^ 


fP/>/* 


fp/V* 


0.14581 


24 3 x 48 


0.4033(62) 


0.350(35) 


0.900(14) 


0.693(70) 


0.14827 


24 3 x 48 


0.3598(65) 


0.313(32) 


0.883(15) 


0.681(69) 


0.15008 


24 3 x 48 


0.3159(79) 


0.274(28) 


0.861(17) 


0.664(68) 


0.15096 


24 3 x 48 


0.291(10) 


0.253(27) 


0.844(21) 


0.650(67) 


0.151959 


24 3 x 48 


0.275(20) 


0.239(30) 


0.845(36) 


0.651(71) 


0.152496 


24 3 x 48 


0.24(16) 


0.21(14) 


0.74(11) 


0.57(10) 



Table 25. Decay constants for the SU(2) theory. 
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K 


iV 3 x N t 


K at /V^ 


F v /yft 


fl at /V^ 


fp/V* 


0.15002 


24 3 x 48 


0.4132(58) 


0.349(35) 


0.9402(83) 


0.686(69) 


0.1522 


24 3 x 48 


0.3727(51) 


0.315(32) 


0.9143(85) 


0.667(67) 


0.1538 


24 3 x 48 


0.3325(64) 


0.281(29) 


0.8835(91) 


0.644(65) 


0.15458 


24 3 x 48 


0.3083(78) 


0.261(27) 


0.865(10) 


0.631(64) 


0.155638 


24 3 x 48 


0.2742(78) 


0.232(24) 


0.849(18) 


0.619(63) 


0.15613 


24 3 x 48 


0.246(14) 


0.208(24) 


0.893(55) 


0.652(77) 



Table 26. Decay constants for the SU(3) theory. 



K 


N 3 x N t 


K at /V° 


iVv^ 


'iflyfi 


fp/V* 


0.15184 


24 3 x 48 


0.4305(37) 


0.360(36) 


0.9807(64) 


0.701(70) 


0.154 


24 3 x 48 


0.3857(39) 


0.323(32) 


0.9526(67) 


0.681(68) 


0.15559 


24 3 x 48 


0.3417(46) 


0.286(29) 


0.9206(77) 


0.658(66) 


0.15635 


24 3 x 48 


0.3227(34) 


0.270(27) 


0.8998(92) 


0.643(65) 


0.157173 


24 3 x 48 


0.2934(49) 


0.245(25) 


0.870(14) 


0.622(63) 


0.15764 


24 3 x 48 


0.2710(71) 


0.227(23) 


0.855(25) 


0.611(64) 



Table 27. Decay constants for the SU(4) theory. 



K 


7V S 3 x N t 


F^/y/o 


F v /yfc 


'iflyfi 


fp/V* 


0.15205 


24 3 x 48 


0.4350(31) 


0.363(36) 


0.9702(50) 


0.688(69) 


0.15426 


24 3 x 48 


0.3909(32) 


0.326(33) 


0.9372(52) 


0.665(67) 


0.15592 


24 3 x 48 


0.3484(37) 


0.290(29) 


0.8984(57) 


0.637(64) 


0.15658 


24 3 x 48 


0.3291(42) 


0.274(28) 


0.8788(65) 


0.623(62) 


0.157548 


24 3 x 48 


0.2951(39) 


0.246(25) 


0.845(10) 


0.600(60) 


0.158355 


24 3 x 48 


0.269(14) 


0.224(25) 


0.815(55) 


0.578(70) 



Table 28. Decay constants for the SU(5) theory. 



K 


7V S 3 x N t 


F™-/yft 


F v /yfc 


'iflyfi 


fp/V* 


0.15264 


24 3 x 48 


0.4374(25) 


0.364(36) 


0.9903(36) 


0.698(70) 


0.15479 


24 3 x 48 


0.3944(26) 


0.328(33) 


0.9600(35) 


0.677(68) 


0.15636 


24 3 x 48 


0.3536(30) 


0.294(30) 


0.9248(37) 


0.652(65) 


0.15712 


24 3 x 48 


0.3302(35) 


0.275(28) 


0.9031(42) 


0.637(64) 


0.158051 


24 3 x 48 


0.2987(34) 


0.248(25) 


0.8458(95) 


0.597(60) 


0.158845 


24 3 x 48 


0.2676(91) 


0.222(24) 


0.776(36) 


0.548(60) 



Table 29. Decay constants for the SU(6) theory. 
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K 


iV 3 x N t 


K at /V^ 


F*/yft 


fP/>/* 


fp/V* 


0.152816 


24 3 x 48 


0.4386(23) 


0.364(36) 


0.9918(43) 


0.697(70) 


0.154967 


24 3 x 48 


0.3949(24) 


0.328(33) 


0.9592(44) 


0.675(68) 


0.156547 


24 3 x 48 


0.3529(27) 


0.293(29) 


0.9214(48) 


0.648(65) 


0.157339 


24 3 x 48 


0.3273(31) 


0.272(27) 


0.8961(55) 


0.630(63) 


0.158273 


24 3 x 48 


0.2906(47) 


0.241(24) 


0.8602(84) 


0.605(61) 


0.159062 


24 3 x 48 


0.249(16) 


0.206(25) 


0.814(28) 


0.573(60) 



Table 30. Decay constants for the SU(7) gauge theory. 



K 


N 3 x N t 


F l f/^ 


F v /yfc 


'iflyfc 


fp/V* 


0.15298 


12 3 x 24 


0.4456(31) 


0.369(37) 


1.0055(92) 


0.704(71) 


0.15521 


12 3 x 24 


0.4045(35) 


0.335(34) 


0.986(10) 


0.690(69) 


0.15684 


12 3 x 24 


0.3610(44) 


0.299(30) 


0.958(11) 


0.670(68) 


0.15755 


12 3 x 24 


0.3379(53) 


0.280(28) 


0.939(13) 


0.657(66) 


0.158531 


12 3 x 24 


0.307(12) 


0.254(27) 


0.902(17) 


0.631(64) 


0.15931 


12 3 x 24 


0.296(28) 


0.245(34) 


0.851(48) 


0.595(68) 



Table 31. Decay constants for the SU(17) theory. 
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